Finite-size effect and the components of multifractality in financial volatility 
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Many financial variables are found to exhibit multifractal nature, which is usually attributed 
to the influence of temporal correlations and fat-tailedness in the probability distribution (PDF). 
Based on the partition function approach of multifractal analysis, we show that there is a marked 
finite-size effect in the detection of multifractality, and the effective multifractality is the apparent 
multifractality after removing the finite-size effect. We find that the effective multifractality can 
be further decomposed into two components, the PDF component and the nonlinearity component. 
Referring to the normal distribution, we can determine the PDF component by comparing the 
effective multifractality of the original time series and the surrogate data that have a normal distri- 
bution and keep the same linear and nonlinear correlations as the original data. We demonstrate 
our method by taking the daily volatility data of Dow Jones Industrial Average from 26 May 1896 
to 27 April 2007 as an example. Extensive numerical experiments show that a time series exhibits 
effective multifractality only if it possesses nonlinearity and the PDF has impact on the effective 
multifractality only when the time series possesses nonlinearity. Our method can also be applied 
to judge the presence of multifractality and determine its components of multifractal time series in 
other complex systems. 

PACS numbers: 05.45.Df, 05.45.Tp, 89.65.Gh 



I. INTRODUCTION 

Turbulent flows exhibit multifractal nature. Two main 
pictures of multifractal analysis in the turbulence litera- 
ture deal with the velocity fluctuations using the struc- 
ture function approach and the energydissipation rates 
using the partition function approach [l| . The difference 
between the two pictures is remarkable and there are a 
lot of efforts trying to relate the two sequences of scaling 
exponents. In addition, the power-law dependence of the 
structure function is more or less not perfect. 

It is well known that there are many similarities be- 
tween turbulent flows (velocity time series) and financial 
markets (equity prices) @, 0, 0, 0]. For instance, mul- 
tifractality is one of the most elusive stylized facts in 
financial markets @, 0, 0| . Many different methods have 
been applied to characterize the hidden multifractal be- 
havior of different financial variables, such as the fluc- 

M, the structure function 
ilSlili, the 
multifractal detrended fluctuation analysis (MF-DFA) 
[HI, H3, [H, [13j2lj2^j23 J2^, J2^, ]3^| , the partition func- 
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tion method 
the multiplier method 
approaches [H, H l45l, |4 

from the multifractal analysis of turbulence data 



M |35|, |36|, |37|, |3j, |39|, |4fl |41| , 

2j, and the wavelet transform 
I, some of which are borrowed 

When 
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the return time series is concerned, the structure function 
method and the MF-DFA method apply. For volatility 
time series, the partition function approach and the mul- 
tiplier approach can be utilized. 

The extracted multifractal nature has potential appli- 
cations in financial engineering. Some researchers report 
that the width of the estimated multifractal spectrum is 
correlated to the price fluctuation in the future and thus 
can be used to predict price fluctuations [H, [H, [4l[ . In a 
similar vein, the so-called multifractal volatility has been 
introduced to measure stock market risks, which can be 
used to estimate the dynamic Value-at-Risk of an asset 
[38l | . There is also empirical evidence suggesting that 
multifractal analysis can be used to quantify the degree 
of inefficiency of markets in the sense that more devel- 
oped stock markets have weaker multifractality [I?], HH 
and an emerging market evolves to be more efficient with 
narrowing singularity width [49| . 

However, understanding the origin of multifractality 
in financial markets and its components is still a subtle 
and open problem. For instance, standard multifractal 
analysis of an exactly monofractal financial model gives 
artificial multifractal behaviors (5(j. It is usually argued 
that the fat tails and the long-term power-law correla- 
tions are two possible sources of multifractal nature in fi- 
nancial time series plj . It is obvious that possessing only 
linear correlations is not sufficient for the presence of mul- 
tifractality and a nonlinear process with long-memory is 
necessary to have multifractality [5l| . Numerical inves- 
tigations show that the reported multifractal nature in 
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many real cases is stemmed from the large fluctuations 
of prices ■ 

Empirically, the problem has been studied based on the 
multifractal detrended analysis of financial returns [53| . 
One conventional measure for quantifying the degree of 
multifractality is the width of singularity spectrum 

Aa = a max - a min . (1) 

Alternative measures are also used in some cases [U H3] . 
In this way, one can quantitatively determine the contri- 
bution of the correlation and PDF components [Hi], • 
In order to investigate the influence of temporal cor- 
relations in return series, one can randomly shuffle 
the orig inal data and compare their singularity spectra 
HI HIT H HI H!> EH- All these studies show that 
the shuffled data have non-shrinking singularity width 
AasF, comparable to the original width Aa. These ob- 
servations imply that the heavy-tailed distribution of the 
returns has a crucial impact on the singularity width. To 
understand the impact of the distribution, one can ei- 
ther remove large returns [57| or generate surrogate data 
having a Gaussian distribution while keeping the linear 
correlation of the original data [1^, [3(| EH ■ In addition, 
the temporal correlations can be further separated into 
linear and nonlinear components, which have different 
effects on the multifractal spectrum |53[ . 

In this work, based on the partition function method 
of multifractal analysis, we propose that the apparent 
multifractality can be decomposed into three components 
caused by the nonlinear correlation, the linear correlation 
(long-term power-law memory) and the fat-tailed PDF, 
which can be characterized by the following expression: 

Aa = Aa NL + Aa L M + Aa PD F- (2) 

We find for the first time that the linear correlation 
component Aolm is the outcome of the finite-size effect 
and thus the intrinsic multifractal nature is character- 
ized by the effective multifractality Aa G ff composing of 
the nonlinearity component AaNL and the PDF compo- 
nent AapDF- These properties are very different from 
those of the MF-DFA results in Ref. [Ij. 

We will develop a systemic procedure to quantitatively 
determine the components of multifractality. The anal- 
ysis is carried out based on the daily volatility data of 
the Dow Jones Industrial Average (DJIA) from 26 May 
1896 to 27 April 2007 (totally 30147 trading days). The 
logarithmic returns is defined as 

r(t)=InP t -lnP t _i, (3) 

where P t is the price at time t, and the volatility is com- 
puted as the absolute value of r(t): 

v(t) = \r(t)\. (4) 

This simple definition is widely adopted in the Econo- 
physics community. 



This paper is organized as follows. Section HT1 gives a 
brief description of the partition function approach for 
multifractal analysis and its numerical implementation. 
Section [Ml designs a numerical procedure to confirm the 
presence of finite-size effect in the detection of multifrac- 
tality in time series with a prescribed probability distri- 
bution. In Sec. IIVI the three multifractal components 
are determined for the daily DJIA volatility, where the 
nonlinear correlation component and the fat-tail compo- 
nent are determined by comparing with surrogate time 
series with normal distributions. The impact of the fat- 
tailcdness of the probability distribution on the apparent 
multifractality is further investigated in Sec.[V]using Stu- 
dent's t distributions and the Wcibull distributions. We 
find that the fat-tailed PDF per se has no effect on the 
multifractality and its impact is due to its coupling with 
the nonlinear correlations. We finally provide some con- 
cluding remarks in Sec. IVII 



II. MULTIFRACTAL ANALYSIS BASED ON 
PARTITION FUNCTIONS 

The multifractal analysis in this paper is based on 
the partition function formalism [59l |. Our procedure is 
slightly different from the conventional algorithm in the 
calculation of the partition functions. Consider a time se- 
ries {vi : i = 1, 2, • • • , N}. For a given integer time scale 
£, we choose randomly m integers {j'j : i = 1, 2, • ■ ■ , m} 
uniformly distributed in [1, N — £+ 1]. Then the measure 
enclosed in the interval +£ — 1] is calculated 

£-1 / N 

Ei(£) = J2 v n+k /!>*, (5) 

fc=0 / fe=l 

The q-th partition function is calculated as follows 

N - 

M S) = —^ E ^- (8) 

i— 1 

If the time series is power-law long-term correlated, we 
have 

M q {£) ~ {£/NY (q) . (7) 

If the scaling exponent r(q) is a nonlinear function of q, 
the time series possess multifractal nature. The singular- 
ity strength and its spectrum can be computed according 
to the Legendre transform of r(q) [f39[, 

a(q) = r'{q) and f(a) = qa — t. (8) 

The advantage of using Eq. ([6]) is twofold. First, it gives 
better statistics when £ is close to N. Second, it does 
not need to deal with the situation in which N/t is not 
integer. 
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III. FINITE-SIZE EFFECT 

The finite-size effect in the detection of multifractal- 
ity has been documented for several uncorrelated time 
series. One example is time series having an exponen- 
tial distribution, whose <?-order moment can be derived 

as m 



M q (l) =T(£ + q)L/r{£+l). 



(9) 



When the length L of the time series is small, a spurious 
multifractal spectrum is detected, where the singularity 
width Aq is significantly larger than zero. Only for suf- 
ficiently long time series with L ^ and for I ^ q, the 
moment is approximated by M q (£) ~ It follows 

immediately that t = q — 1, a = 1 and f(a) = 1. There- 
fore, the finite-size effect results in a spurious detection 
of multifractality from a monofractal signal if the length 
of the signal is short or the "scaling range" locates at 
small scales. 

Another example is exponentially truncated Levy 
flights with characteristic parameter 7, which exhibit 
a bifractal behavior such that r(q) = g/7 — 1 when 
< q < 7 and r(q) = when (7^7 [61(. The multifractal 
spectrum f(a) shrinks into two points (0, 0) and (I/7, 1). 
Extensive numerical experiments using uncorrelated time 
series obeying q-Gaussian distributions with different tail 
exponents unveil a convergence to monofractalilty in the 
Gaussian attraction basin and to bifractality in the Levy 
attraction basin [631, which is consistent with the ana- 



lytic results for exponentially truncated Levy flights [61 



More generally, according to the law of large numbers, 
the sum of I i.i.d. random variables with mean \x and 
standard deviation a converges to ly, with much smaller 
fluctuations [63[ . Let u, > be a sequence of independent 
random variables with mean \i and variance a 2 and 



Denoting that 



and 



we have 



Xi = Vi — n 



X = Y^Xi, 



V = X + lfj,. 



(10) 



(11) 



(12) 



(13) 



According to the law of large numbers [63j , X converges 
to 0, if 

/riPW-O-d/^Xoo, (14) 



where F{x) is the cumulative distribution of x. Hence, 
V converges to In addition, according to the Central 
Limit Theorem, the distribution of X is Gaussian: 



P(X) 



1 



■ exp 



(X - in? 



2£<j 2 



(15) 



We have 



M q (t) = jj V<*P(V)dV 



N 



{x + inYP{x)dx 



NP- 1 J P(X)dX 



(16) 



where P(X)dX = P(V)dV and X < Ip for most X are 
used. It follows immediately that t = q — 1, a = 1 and 
f(a) = 1. It indicates that the time scries is monofractal 
and its "singularity spectrum" shrinks to a single point 
(a, f) — (1, 1). A direct consequence is that any degree of 
multifractality observed in real uncorrelated time series 
is caused by a finite-size effect. 

However, the impact of linear long-term correlations 
on the finite-size effect has not been investigated. It is 
well-known that financial volatility exhibits strong long- 
term memory @, 113. The Hurst index of the daily DJIA 
volatility is H = 0.80. In order to show that there is a 
finite-size effect in the detection of multifractality intro- 
duced by linear correlations in the volatility time series, 
we generate surrogate data that have the same proba- 
bility distribution as the original time series of volatility, 
which can be done with the transformation method [65| . 
We first construct the empirical cumulative distribution 
F(v) of volatility, which is the occurrence frequency of 
volatility v(t) less than v. A sequence of random num- 
bers {xi : i = 1,2, •■■ ,L} arc drawn from a uniform 
distribution. In order not to introduce very large val- 
ues in the surrogate data of volatility caused by the edge 
effect, we perform the following map 



Vi 



(xi — min{xj})(l — min{_F}) 



max{xj} — minlxj} 



min{^}. (17) 



The sequence of numbers {zi 



1,2, ■• ■ , L} are calcu- 



lated according to the following transformation 

Zi = F-\y l ), (18) 

which can be done numerically by linear or spline inter- 
polations. The sample {zi : i = 1, 2, • ■ • , L} has the same 
distribution as the volatility {v(t) : t = 1, 2, • ■ ■ , N}. 

We then introduce long memory (linear correlations) 
in the time series {z{\ using an improved am plit ude ad- 
justed Fourier transform (IAAFT) algorithm [66j, which 
is based on a simple iteration scheme and an improved 
version of the amplitude adjusted Fourier transform al- 
gorithm [13. The volatility data {v(t) : t = 1, 2, ■ • • , N] 
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are sorted resulting in a new sequence {sat}, and we ob- 
tain the squared amplitudes of the Fourier transform of 

{sjv}, denoted as {S%}. The initial sequence {s^ '} of 
the iteration is a random shuffle of {s^r}. In the i-th 

iteration, the squared amplitudes {Si } of the Fourier 

transform of {s}^ } are obtained and replaced by {5^}, 
which are transformed back, and then the the resulting 
series are replaced by {sjv} but keeping the rank order. 

In our numerical experiments, we have investigated 
different Hurst indexes with ff = 0.1,0.2,---,0.9. The 
length L of the surrogate data ranges from 10 3 to 10 7 . 
Figure Q] shows the dependence of M g (f) 1 ^ <?_1 ^ as a func- 
tion of £ for different surrogate time series with different 
linear correlations. According to Eq. ([7]), the slope of lin- 
ear lines in Fig. [1] is r(q)/(q — 1). For monofractal time 
series, these lines should be parallel and have identical 
slopes. For short time series, the scaling range locates on 
the left of plot and the lines are not parallel. Only when 
the time series is sufficiently long, the monofractality can 
be detected. More data points are needed to reach a de- 
cisive conclusion for strongly correlated time series with 
large Hurst index H . 



which accounts for the contribution of the nonlinear cor- 
relations and the fat tails to the apparent multifractal 
spectrum width. 




FIG. 2: (Color online) Quantification of the finite-size effect 
in the detection of multifractality. 




1/(9-1) 



func- 



FIG. 1: (Color online) Dependence of M q (£) 
tion of £ for different surrogate time series with different linear 
correlations. The Hurst index ranges from 0.9 to 0.1 from the 
top to the bottom. 



For each surrogate time series with Hurst index H and 
size L, a multifractal spectrum is determined within the 
scaling range [L/60, L/3]. Figure [2] illustrates the finite- 
size effect in the detection of multifractality. For a given 
H, the singularity width Act is a decreasing function of 
the series length L. It is shown that the observed non- 
vanishing singularity spectrum in the surrogate data with 
only linear correlations is just the outcome of a finite-size 
effect, 

Aa L M = Aa F SE- (19) 
Therefore we can define an effective singularity width 

Acteff = Acv - Aa LM = Aa NL + Aa PDF (20) 



According to Fig. [2l there is a power-law dependence 
between the singularity width Aa(H, L) and the length 
L for each Hurst index H when L is not large: 



Aa(H, L) ~ L a 



(21) 



where a is a function of H. For each H, we calculate a in 
the scaling range [1000, 138950]. The dependence of the 
exponent a on the Hurst index H is shown in Fig. [31 We 
find a nice linear relation 



2H-2 



(22) 



except for H = 0.1. It is interesting to note that, for 
H = 0.5, we have a = — 1 and thus Aa(H, L) oc L _1 . 




-1.5 



FIG. 3: Dependence of the power-law exponent a on the Hurst 
index H. 
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We assume that Aa(H, L) can be factorized as follows 

Aa(H, L) = g(H, L)L a . (23) 

Figure [4] plots g(H,L) as a function of L for each H. 
We find that g(H, L) is almost independent of L, in the 
scaling range. For the special case with the Hurst index 
H = 0.5, it is observed that the g(H, L) curve is horizon- 
tal in the whole range of L. 




L 

FIG. 4: (Color online) Dependence of g(H,L) as a function 
of L for different Hurst index H. 

According to Fig. 0J we can assume that g{H, L) is 
independent of L in the scaling range. We calculate the 
average (g(H, L)) of g(H, L) in the scaling range 

(g(H,L)) = - J2 9(H,L) (24) 

Le[1000, 138950] 

where n is the number of scales L in the scaling range. 
Figure[5]plots h\{g{H, L)) as a function of H. We observe 
a linear relation 

ln(g(H, L)) = 10 - 10H, (25) 

which is obtained based on the linear least-squares re- 
gression. 

Combining Eqs. (E2) and ([25]), we obtain 

Aa^^^i-^-^e 10 ^"^ (26) 

when L is not large, ft is interesting to note that 
Aa(H, L) oc L^ 1 for H = 0.5 in the whole range of L. 
Note that this expression (f2"6")l is valid only for not large 
L. Also, in the quantitative assessment of finite-size effect 
in correlated signals, one should use Fig. O rather than 
the above expression. For other applications, one need 
to conduct numerical simulations to produce a finite-size 
effect figure as Fig. [5] based on the sample distribution 
under investigation. 




FIG. 5: Exponential dependence of (g(H,L)) with respect to 
the Hurst index H. 



IV. DETERMINING THE SINGULARITY 
WIDTH COMPONENTS 

We now try to determine the components of the ap- 
parent singularity width Aa. Figure [6] shows the mul- 
tifractal spectrum of the original daily DJIA volatility 
time series. The apparent singularity width is found to 
be Aa = 0.367. 




0.97 0.99 1.01 1.03 1.05 

a 



8.7 0.8 0.9 1 1.1 1.2 

a 

FIG. 6: (Color online) Determination of the apparent singu- 
larity width Aa, the finite-size effect component Aq:fsb, and 
the effective singularity width Aq c h of the daily DJIA volatil- 
ity. The inset is the amplification of the main plot around 
(a,/) = (1,1). 

In order to determine the finite-size effect component 
caused by the linear long-term correlations in the volatil- 
ity series, we use the iteration scheme [6(| to generate 
surrogate data Ulm which have the same distribution 
and linear long-term correlations as the original volatil- 
ity time series, while any underlying nonlinear correla- 
tions have been eliminated. One hundred surrogate time 
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series wlm have been generated and the average multi- 
fractal spectrum is shown in Fig. [5] The finite-size effect 
component is found to be 



AapsE = Aa LM = 0.219 ± 0.038. 



(27) 



For the shuffled data i>sf, the resulting multifractal spec- 
trum shrinks to a point (a, /) = (1, 1) as expected and 
the singularity spectrum width AasF — 0.007 ± 0.001 
is close to zero, which confirms that the fat-tailedness 
of PDF alone cannot produce any multifractality. This 
finding is very different from that in the multifractal de- 
trended fluctuation analysis of financial returns. It fol- 
lows that the effective width of singularity is 



Ad- 



en 



Aa - AapsE = 0.148 ± 0.038, 



(28) 



which is significantly smaller than the apparent singular- 
ity width Aa. 

In order to determine the PDF component AapDF, 
we need to generate surrogate data in which the linear 
and nonlinear correlations are preserved while the orig- 
inal PDF is replaced by a reference PDF such that the 
PDF of the surrogate data has no impact on the result- 
ing multifractality. It is not irrational to assume that the 
PDF component of Gaussian signals is negligible, that is, 
Aa normj pDF = 0. We generate surrogate data r norm for 
the returns r(t) which are drawn from a normal distribu- 
tion 

1 (-noi-m-M) 2 , N 

p(r„ OT m) = -7=^e ^ , (29) 



2na 

where [i = 0.00021 and a = 0.011 (in units of one trading 
day) are the sample mean and standard deviation of the 
daily DJIA returns r(t). The surrogate volatility is thus 
fnorm.SF = |?"norm|j which has no temporal correlation 
and can be regarded as being shuffled. 

The random numbers w n orm,SF are rearranged to have 
the same rank ordering as v(t) so that the surrogate 
time series v norm (t) preserves the linear and nonlinear 
correlations of the original volatility v(t) [H, [H, [69l |. 
The approach is to replace the raw data by random 
numbers drawn from a prescribed distribution, which 
is described as follows. For a given distribution, we 
generate a sequence of random numbers {xo(t) : t = 
1, 2, ■ ■ • , N}, which are rearranged such that the result- 
ing series {x(t) : t = 1, 2, • • • , TV} has the same rank or- 
dering as the volatility series {v(t) : t = 1,2,- •• , TV}. 
In other words, x(t) should rank n in the sequence 
{x(t) : t = 1,2,- •• , N} if and only if v(t) ranks n in 
the {v(t) : t = 1,2,- •■ ,N} sequence [H, M, HI • 

We generate 100 surrogates v norm (t). For each time se- 
ries i>norm(i), we further destroy its nonlinear correlations 
but keep its linear correlations to generate iWm.LM^) 
based on the iterated amplitude adjusted Fourier trans- 
formation algorithm [66| . The average multifractal spec- 
tra for v nolm (t), i>norm,LM(i)> and i>norm,SF(<) are illus- 
trated in Fig. [7l We find that, in the Gaussian surrogate 
case, the apparent singularity width is 



for D norm (t), its finite-size effect component is 

Aa norm ,FSE = Aa norm xM = 0.118 ± 0.028 (31) 

for W n orm,LM(i), and 

Aa no rm,SF = 0.003 ± 0.001 (32) 

for Wnorm.SF(i)- Therefore, the effective width of singu- 
larity for the Gaussian surrogates is 

Ao-norm.eff = Aa norm - Aa norm . FSE = 0.067 ± 0.028. 

(33) 

Taking into account the assumption that Aa normi pDF = 
0, it follows immediately that the two components of the 
effective multifractality of the original volatility are 



AapDF = Aa e 



Aa 



norm.crl 



0.081 ±0.033 (34) 



and 
Aa N L 



^norm,NL 



= Aa norm , cff = 0.067 ± 0.028. (35) 



We note that, in the above framework, the nonlinearity 
component does not change when the PDF is replaced. 




FIG. 7: (Color online) Determination of the PDF component 
AapDF and the nonlinearity component Aqnl of the original 
volatility with reference to Gaussian surrogates. The inset is 
the amplification of the main plot around (a, f) — (1, 1). 



V. IMPACT OF PDF 

In order to further confirm that the PDF of the time 
series has crucial influence on the multifractality, we in- 
vestigate two families of distributions with fat tails. The 
first one is a family of Student's t distributions 



- 7+1 ' 
* 2 / 



(n - m) 2 

7 



"(7+l)/2 



(36) 



Aa„ 



0.185 ±0.003 



(30) 



which have power-law tails with exponent 7. In our sim- 
ulation, 7 ranges from 3 to 10 with a spacing step of 0.5. 
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For each tail exponent 7, we generate three types of sur- 
rogates for the volatility v(t) denoted as v t (t) = \r t (t)\ 
which has the same temporal correlations as v(t), wt,LM(i) 
which has the same linear correlations as v(t) while any 
nonlinear correlations are eliminated, and i>t,SF(i) which 
is the shuffled data of vt(t). We have generated 100 se- 
ries of Vt(t) which has the same length as v(t). For each 
Vt(t), we generate two surrogates i>t,LM(i) and vt,hM.(t)- 
The multifractal spectrum of each time series is deter- 
mined. 

The apparent singularity width Aa t , the finite-size ef- 
fect component Aa t ,FSE, the effective singularity width 
Aat, e ff, the PDF component Aa t ,PDF are calculated and 
illustrated in Fig. [8] as a function of the tail exponent 7. 
All the quantities decreased with 7, implying that the 
multifractality is stronger when the distribution is fat- 
ter (characterized by smaller 7). We also plot Aa e ff and 
AapDF for the original volatility as horizontal lines, both 
of which intersect with the Aa t ,eff(7) and Ao^pdf^) 
curves at 7 = 3. For 7 = 3, we have Aa t = 0.387± 0.008 
for vt(t), Aa t ,FSE = 0.241 ±0.062 for v t , LM (t), Aa t , S F = 
0.008 ± 0.002 for v i>SF (t), Aa tjCff = 0.145 ± 0.062 and 
Aa t ,pDF = 0.078 ±0.062. This striking feature is consis- 
tent with the inverse cubic law of financial returns [70l |. 
which also holds for daily DJIA returns [7l| . 




FIG. 8: (Color online) Dependence of the singularity width 
Aa as a function of the exponent 7 for surrogate time series 
with Student's t distribution. The error bars are the standard 
deviations for the 100 surrogate series. 



The second one is a family of Weibull distributions 



p(v w u) = Pv^e-"^ 



(37) 



where the shape parameter (3 describes the heaviness of 
the tails and we require that (3 ^ 1. In our simulation, 
(3 varies from 0.4 to 1.0 with a spacing step of 0.1. For 
each /3, 100 time series of the same length of v(t) are 
generated, which are manipulated to generate i> w bi with 
the same rank ordering of v(t) and u w bi,LM only with 
linear correlations. 




3 0.4 0.5 0.6 0.7 0.8 0.9 1 1.1 
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FIG. 9: (Color online) Dependence of the singularity width 
Aa as a function of the exponent (5 for surrogate time series 
with Weibull distribution. The error bars are the standard 
deviations for the 100 surrogate series. 



The results are shown in Fig. [9l We find that Aa w bi, 
Aa wW ,FSE and Aa wW ,sF decrease with /?, while Aa w bi, e ff 
and Aa w bi,PDF are concave functions with respect to f3. 
For small /3 values, Aa w bi, e ff an d Aa w bi,pDF increase 
with (3, which may be caused by the fact that there is 
a lack of statistics so that the fluctuations of Aa w bi,FSE 
are markedly large. Indeed, the error bar decreases with 
increasing (3. It is worth noting that the fluctuations 
of Aat and Aa w bi are much less than the corresponding 
Aat,FSE and Aa w bi,FSE- The large error bars of the PDF 
components and the nonlinearity components are caused 
by the large error bars of the corresponding finite-size ef- 
fect components. In addition, in the case of exponential 
distribution with j3 = 1, we have Aa oxp = 0.318 ± 0.004 
for v cxp (t), Aa cxp ,FSE = 0.198±0.005 for u eX p,LM(t), and 
Aa cxp ,SF = 0.006 ± 0.001 for v cxp , SF {t). It follows that 
Aa e xp!cff = 0.120±0.005 and Aa e x P ,PDF = 0.053±0.005. 



VI. SUMMARY 

We have proposed a method to decompose the appar- 
ent multifractality into three components associated with 
the finite-size effect, the nonlinearity, and the fat-tailed 
PDF. The finite-size effect component can be determined 
by the singularity width of surrogate data with the same 
PDF and linear correlations as the original time series 
after eliminating its nonlinear correlations. The PDF 
component can be calculated as the difference between 
the effective multifractality of the original series and the 
surrogate data drawn from a normal distribution with 
the linear and nonlinear correlations preserved. The im- 
pact of the fat-tailed PDF on the effective multifractal- 
ity is the outcome of the coupling of the PDF and the 
nonlinearity since the shuffled data are monofractal with 
vanishing singularity width. The method proposed in 
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this paper can be applied to investigate the presence of 
multifractality in time series and, if any, determine its 
components. 

It is necessary to clarify the marked differences between 
the current work and Ref. [53|. First, the current work in- 
vestigates financial volatility, while Ref. [53[ investigates 
financial returns, although the same data set of the daily 
DJIA index is used in both papers. Second, the parti- 
tion function approach is used in this work, while the 
multifractal detrended fluctuation analysis is adopted in 
Ref. [53|. Third, the results are quite different. Ref. [53j | 
shows that the multifractal spectrum width AasF of the 
shuffled data is comparable to the apparent singularity 
width Aa, which means that the long-term correlations 
have minor impact on Aa and the fat-tailed PDF plays 
a major role. In contrast, in Figs. and [5] of the 

current paper, we have shown that 

Aa SF « (38) 

for shuffled time series with different PDFs. It is clear 
that 

Aa PDF ± Aa SF - Aa norm , SF , (39) 

where AasF is the singularity width of the shuffled data 
i>sf and Aa normi sF is the shuffled data v n orm,SF- This 



suggests that the contribution of the PDF is coupled 
with the nonlincarity. If there is no nonlincarity, the 
fat-tailedncss of the PDF will not introduce any multi- 
fractality. In addition, with the presence of nonlincarity, 
the PDF does have impact on the effective multifractal- 
ity. 

As a last note, it is well-established that the 
phase-randomized surrogates of heartbeat series from 
healthy subjects have a narrower width of singularity 
(Aa nornlj LM ~ 0.1) compared to its apparent multifrac- 
tal width (Aa ss 0.45) based on the wavelet transform 
method (73 . [73[ , which is in excellent agreement with our 
results. 
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